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Abstract. The problem of calculation of correlation functions in the six- 
vertex model with domain wall boundary conditions is addressed by consider- 
ing a particular nonlocal correlation function, called row configuration proba- 
bility. This correlation function can be used as building block for computing 
various (both local and nonlocal) correlation functions in the model. The 
row configuration probability is calculated using the quantum inverse scat- 
tering method; the final result is given in terms of a multiple integral. The 
connection with the emptiness formation probability, another nonlocal corre- 
lation function which was computed elsewhere using similar methods, is also 
discussed. 



1. Introduction 

One of the most fundamental problems in the theory of integrable models is 
the exact calculation of correlation functions pQ. In recent years, there has been an 
increasing interest (motivated by various mathematical and physical applications) 
in obtaining exact results for correlation functions of statistical mechanics models 
defined on finite lattices and with fixed boundary conditions. Because of the lack 
of translational invariance, the systematic computation of correlation functions for 
these models represents a difficult problem. 

An important example of such a model is the six- vertex model with domain wall 
boundary conditions [2] • The partition function of the model on the finite lattice is 
given exactly in terms of certain determinant [3,4. This formula, known as Izergin- 
Korepin determinant formula, turned out a powerful tool in proving important 
combinatorial results. The current interest in the model is mostly motivated by 
occurrence of the phase separation phenomena (see [5j[6] and references therein). 

Some progress in the calculation of the correlation functions of the six-vertex 
model with domain wall boundary conditions has been achieved when correlations 
are considered near the boundaries [7rU0j. An example of correlation function 
which can be computed away from the boundary is the so-called emptiness forma- 
tion probability Some generalisations have been considered recently in [12] . 
However the problem of a systematic treatment of correlation functions, especially 
when correlations are considered away from the boundaries, is still far from being 
solved. 

To address this problem, in the present paper we introduce a particular non- 
local correlation function, called row configuration probability. This correlation 
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function describes the probability of observing a given configuration of arrows on the 
vertical edges located between two consecutive horizontal lines of the square lattice. 
The row configuration probability can be used as building a block to compute other 
(both local and non-local) correlation functions. In particular, it is closely related 
to the emptiness formation probability. 

The row configuration probability besides being interesting for evaluation of 
other correlation functions, is also interesting on its own right: there are analogue 
in the context of phase separation for dimers models |13j , and of enumerative com- 
binatorics P3]. 

To compute the row configuration probability we consider the inhomogeneous 
version of the model and formulate it in the framework of the quantum inverse 
scattering method (QISM) [TS] (for a survey, see [T]). We use the fact that the 
row configuration probability can be represented as a product of two factors. For 
computing the first factor we use a side result of paper [IB] , while for the second one 
we use the technique developed in papers [H[T0l[Tl]. For the homogeneous model 
both factors are represented in terms of multiple integrals. 

To demonstrate how these results can be used for computing other correlation 
functions, we discuss here the connection with the emptiness formation probability. 
The computation is based on performing certain sums and integrals, and making 
use of identities involving antisymmctrisation of multi- variable functions. 

The paper is organized as follows. In the next section after recalling basic facts 
about the model we set up the considered problem in terms of QISM objects. The 
derivation of the row configuration probability is given in sections 3 and 4. The 
relation with the emptiness formation probability is discussed in section 5. 

2. Six-vertex model, domain wall boundary conditions, and row 
configuration probability 

We consider the six-vertex model on a square lattice formed by intersection of 
N horizontal and N vertical lines (an N x N lattice), with special fixed bound- 
ary conditions called domain wall boundary conditions. Recall that the six-vertex 
model is a model in which local states are arrows pointing along edges of the lat- 
tices; the allowed arrow configurations are subject to the 'ice-rule': each vertex 
should have the same number of incoming and outgoing arrows. The Boltzmann 
weights are assigned to the six possible vertex configurations of arrows allowed by 
the ice-rule, and in the model invariant under reversal of all arrows there are three 
different Boltzmann weights, usually denoted a, b, and c. The domain wall bound- 
ary conditions mean that all arrows on the left and right boundaries are outgoing 
while all arrows on the top and bottom boundaries are incoming, see figure [TJ 

To use QISM in calculations we consider the inhomogeneous version of the 
model, in which the weights of the vertex being at the intersection of ath vertical 
line (enumerated from the right) and kth horizontal line (enumerated from the top) 
are a a , k = a(\ a ,v k ), b a , k = b(\ a ,v k ), and c a , k = c, where 

a(A, v) = sin(A — v + rj), 6(A, v) = sin(A — v — rf), c = sin2ry. (2.1) 

The parameters Ai, . . . , Xn are assumed to be all different; the same is assumed 
about vi, ... , i/jv- The parameter 

n 2 -I- 7) 2 — c 2 

A = — ■ — = cos2?7 (2.2) 
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Figure 1. The six vertices and their weights (left), and N-by-N 
square lattice with domain wall boundary conditions (right). 



takes the same value for all vertices, that ensures integrability [17] . The partition 
function is defined as follows 

N 

Zn = J2 II W °AC), (2.3) 

C a,k=l 

where w a ,k(C) takes values w Q ,fe(C) = a a k, b a k, c a k, depending on the configuration 
C. Clearly, Zn — Zn(\\, . . . , Ajv; fi, . . . , vn) where Ai, . . . , Ajv and v\, . . . , vn can 
be regarded as 'variables'; parameter rj has the meaning of a 'coupling constant' 
and it is often omitted in the notations. After QISM calculations, the homogeneous 
model quantities (e.g., partition function) can be obtained from the inhomogeneous 
ones upon setting X a = A (a = 1, . . . , N) and Vk = v (k = 1, . . . , N), where, with 
no loss of generality, one can further put v — 0, see (|2.ip . We shall refer to this 
procedure as homogeneous limit. 

We now define the main objects of QISM in relation to the model. First, let us 
consider vector space C 2 and denote its basis vectors as the spin-up and spin-down 
states 

It) ^ (J) , II) = (?) ■ (2.4) 

To each horizontal and vertical line of the lattice we associate vector space C 2 . We 
also use the convention that upward and right arrows correspond to the 'spin up' 
state while downward and left arrows correspond to the 'spin down' state. 

Next, to each vertex being intersection of the crth vertical line and the /cth 
horizontal line we associate the operator L a ^{\ a , Vk) which acts nontrivially in the 
direct product of two vector spaces C 2 : in the 'horizontal' space T-Lk = C 2 (associ- 
ated with the fcth horizontal line) and in the 'vertical' space V a = C 2 (associated 
with the ath vertical line). Referring to the scattering matrix picture, the arrow 
states on the top and right edges of the vertex can be regarded as 'in' indices of 
the L-operator while those on the bottom and left edges as 'out' ones, that gives 

L a , k {\ a ,v k ) = a Q , fc ^-^ + b a . k ^-^ + c a>k {T-<rt + T a<^k)- ( 2 -5) 

Here T l a and a\ (I = +, — , z) denote operators acting as Pauli matrices in V a and 
Hk, respectively, and identically elsewhere. 
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An ordered product of L-operators along a vertical (or horizontal) line of the 
lattice corresponds in QISM to a monodromy matrix. To construct, e.g., the vertical 
line monodromy matrix, it is useful to think of the L-operator as a 2 x 2 matrix 
acting in space V Q , with the operator entries acting in the space T-L\ : ....n — ®k=i W-k, 
i.e., 



i a ^ v )—JL +b (A,i/)_^-t c-a h 



(2-6) 

Here the subscript indicates that this is a matrix in V a . The ordered product along 
the crth vertical line is the 'vertical' monodromy matrix: 

= fAl... <N (\ a ) Bl N (X a ) 



(2.7) 

[VJ 



The operators n(X) = ^i' jv(^> • • • > ^jv), etc, act in "Hi,. ...at and they are 
independent of a. Each of these operators corresponds to a vertical line of the 
lattice, with the top and bottom vertical arrows fixed. 

Similarly, one can consider the 'horizontal' monodromy matrices, 

Tki^k) = L Ntk (\ N ,v k ) ■ ■ ■ L 2 ,k(X2,Vk)Li tk {\i 1 Vk) 
= Ml...,jvK) B?_ N (v k )\ 

where operators N (v) — A^ N {v; Ai, . . . , \n)i etc, act in Vi,...,jv '■— ®a=i Hx- 
Each of these operators correspond to a horizontal line of the lattice, with the 
rightmost and leftmost horizontal arrows fixed. 

The importance of the monodromy matrix operator entries is that they obey a 
quadratic algebra, called the algebra of monodromy matrix or Yang-Baxter algebra 
PQ. The algebra involves in total 16 commutation relations, and in the following 
we will need some of these commutation relations, namely 

B{X) -B(A') = -B(A') B(X), C{X)C(X') = C(X')C(X), (2.9) 

and 

A(X) B(X') = /(A, A') B(X') A(X) + g(X', A) B(X) A(X'), (2.10) 
where A(X) = A^[ jy(A), etc, and functions /(A', A) and g(X',X) are 

tl\> \\ sin(A - A' + 2rj) , , sm2r] 

f{X , A = — — — , g A , A = — — . (2.11 

sm(A — A') sm(A — A') 

Exactly the same relations are also valid for A(y) = A^ etc (after replacing 

A i->- v and A' i-> v 1 in ()2.9[) and (|2.10jl ). For the full list of commutation relations, 
as well as for details of their derivation, we refer to chapter VIII of book [1] . 

Now we are ready to formulate the model in the framework of QISM. Denoting 
by \tk) and the basis vectors of space Hk, let us introduce states 

IC"> : = ® It*>. IC,JV>:= ® |+*>- (2-12) 

k—1 k—1 



These states are the 'all spins up' and 'all spins down' states in the space Hi at, 

respectively. Taking into account that the domain wall boundary conditions select 
for ath vertical line the operator BY at(A q ), we can write the partition function 
as the matrix element: 

N 

Z N (X U . ..,X N ;v 1 ,...,v N ) = (^l..., N \ J] B L~A X <*) IC,iv) • (2-13) 

ct=l 

We also recall that N (X) = B^ n {X;v\,...,vn). 

Essentially in the same way, one can construct the partition function consider- 
ing operators associated with the horizontal lines. Denoting by |f^) and \\%) the 
basis vectors of V a , we can introduce states 

|^...,iv) := ® |t*>, := I |4£>, (2.14) 

a— 1 a— 1 

which are the 'all spins up' and 'all spins down' states of space Vi....,n- The partition 
function reads: 

N 

Z N (Xi,...,X N ;u 1 ,... t u N ) = <C..,Jv|II C 'b...,iv( I/ fc)|C..^>, (2-15) 

fc=i 

and we recall that Cf- N (v) = Cf N (v; Ai, . . . , Ajv)- 

The partition function is known to be given by Izergin-Korepin determinant 
formula (see [2H1]) 

7 na = inf = i a (Aq,^)fe(A a ^fc) 

where d(X, A') := sin(A — A') and 

y(A, v) — C (2.17) 

a[X, v)byX, v) 

while a(X, u), b(X, v) and c are defined in (|2.ip . For the original proof of (|2.16l) sec 
[I] ; an alternative derivation of this formula can be found in [HUH] ■ 

In the homogenous limit, i.e., when Ai = • • • = A^v = A and v\ = • • • = i>n = 0, 
expression ()2.16j) becomes 

Zn{\ . . . , A;0, . . . ,0) = [sMX jT ( { "t V)]N2 K dat {flTMA)} (2-18) 

Hn=l ( n! ) K«,fe<iV 

where </?(A) := y(A, 0). Below we often use simplified notations for the homogeneous 
model quantities, e.g., writing Zpj for Zn(X, . . . , A; 0, . . . , 0), and so on. 

Let us now turn to the row configuration probability. To define this quantity it 
is useful to mention first that in the six-vertex model with domain wall boundary 
conditions all configurations are such that on the sth row (i.e., on the N vertical 
edges between the sth and the (s + l)th horizontal lines, counted from the top, 
in our conventions) there are exactly s arrows pointing up. It is therefore natural 
to study the probability of observing a given configuration of arrows on a given 
row, or row configuration probability, for short. Namely, we denote by H^ 1 ^'"' 7 "^ 
the probability that the s up-arrows of the sth row are exactly at the positions 
t*i, . . . ,r s (counted from the right), see figure [5] 

Since the row configuration probability describes generic configurations of the 
model, it can be used as a building block to compute other correlation functions. 
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Figure 2. Row configuration with s up arrows between sth and 
(s + l)th horizontal lines from the top, and the resulting partition 
functions for the top and bottom portions of the original lattice. 



For example, by properly summing over positions of up-arrows, one can recover the 
so-called the emptiness formation probability studied in This connection is 

discussed in section 5. 

To compute the row configuration probability, we separate the original N x N 
lattice into two smaller lattices: an upper lattice, with s horizontal and N vertical 
lines, and a lower lattice, with N—s horizontal and N vertical lines. We shall denote 
Z*° p rs and Z^° l Ts the partition functions of the six-vertex model on the upper 
and lower sublattices, respectively (see figure[2]). The row configuration probability 
is essentially given as a product of the partition functions of the six- vertex model 
on these two smaller lattices, 

ytop yhot 

jjfr r.) = n.-.r ri ,..., r , _ (2 19) 

Our main goal in the present paper is therefore the derivation of some useful repre- 
sentations for the partition functions Z*° p and Z^°* ra . Specifically, we provide 
multiple (s-fold) integral representations for these quantities. 

In terms of QISM objects, the partition functions on the upper, sxN, sublatticc 
can be written similarly to representation (|2.15|) . as follows: 

s 

KZ..r = (^l..,N\T- ■ ■ ■ T-T- J] Cf_M |^... iJV ) , (2.20) 

fe=l 

where tJ {j = 1, . . . , N), as above, denote Pauli matrices acting in spaces Vj. 

To write the partition function of the lower, (N — s) x N , sublattice as a matrix 
element, let us define vectors: 

!&,...,*>= ® |t£>, ld,...,*> = ® \ik)- ( 2 - 21 ) 

fc=s+l k=s+l 
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These are 'all spins up' and 'all spin down' states of the space T-L s +i,...,n ■ Cor- 
respondingly, let us consider the matrix elements of the 'truncated' vertical mon- 
odromy matrix given as the product L Q! jv(A a , vn) • • ■ -^a,s+i(A Q , v s +i)- These ma- 
trix elements are operators jy(A a ) = N (X a ; 24+1, . . . , fjv), etc, acting 
in H s +i,...,n ■ The partition function Z^°* r< can be written as 

N r,-l 

^.-,r. = yC+Wv| II S(A a )-A(A r J. J] B[X a ) 

a— r s +1 a— r s _ i + l 

r 2 -l ri— 1 

x---xA(X r2 ) [] S(A a )-A(A ri )- J] B(\ a )\ftY +h ..., N ), (2.22) 

a— ri + 1 a— 1 

where A(X) := ^(A) and B(X) := B S V +1 ,...^(A). 

3. Calculation of ^ 

The matrix element in (|2.20j) can be formally evaluated (as a function of 
Ai, . . . , Ajv and v\, . . . , v s ) using the equivalence of the algebraic and coordinate 
versions of the Bcthc Ansatz. This equivalence was first explicitly proved, as a side 
result, in [16] (see appendix D of that paper); see also book p], Chapter VII. 

For simplicity, we start directly from the case where parameters Ai, . . . , Xn are 
already taken to the same value A, but the parameters z/i, . . . , v s are left arbitrary 
(and not equal to each other). Equation (D.4) of reference [16] in such a case implies 

iJV-1 



where 



x^Hfnc n (tp i t Pt ,-2At Pj +i), (3.i) 

PgSl s j=l l^j<k^s 

t k := b ^p\, (3.2) 
a(A, v k ) 

and the sum is taken over elements of the symmetric group Q s , i.e., permutations 
P ; 1, . . . , s \-¥ Pi, . . . , P s , with [P] denoting the parity of P. Clearly, the expression 
standing in the second line in (13.11) is exactly the s-particle coordinate Bethe Ansatz 
trial wave-function. 

To study the homogeneous limit of (|3.1[) in the remaining set of parameters, we 
first transform slightly this expression. Let us set tk — t + x k (k = 1, . . . , s) where t 
is an arbitrary parameter and the new parameters xi, . . . , x s arc all different. Using 
the fact that for a function /(x), regular near point x = t, one can always write 
f(t + x) = exp(xd z )f(t + z)\ z= o, wc bring (|3.1I) to the form 

Z r°l..,r s = cS f[ KA^fc)]^ 1 II ^h^ 1< f k t <A e ^^ d ^ 

k—1 l^j<k^s 3 ^ ' ^ 

s 

x Ylit + zjp- 1 Yl [(t + z j ){t + z k )-2A(t + z j ) + l] 

3=1 l<i<fe<s zx=—=z s =Q 

(3.3) 

which simply represents an equivalent way to write 
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Let us now consider the homogeneous limit in the parameters v\, . . . , v s . Since 
t is arbitrary, we can perform this limit such that tx = • • • = tfe = t in the limit, 
and put t — b/a, where a and b are the homogeneous model weights, see (|2.1I) . We 
thus have to consider (|3.1[) at x\ = ■ ■ ■ = x s =0. The limit in (|3.3[) can be done 
using the relation 



Yll^j<k^s( X k x j) 



det ( /■ \^ 



(3.4) 



a:i — a; s — 

Reexpressing the values of derivatives at Z\ = ■ ■ ■ = z 8 = as residues, we obtain a 
multiple integral representation 



ytop 

ri,.. 



Co Co 



3=1 



["I [(t + «i)(t + *fc)-2A(i + «,•) + !] 



d s z 
(2i7r) f 



(3.5) 

7 177' I ' 

X<i<A<« 

Here Co is a small, simple, closed, positively-oriented contour enclosing point z = 0. 
Evaluating the Vandermonde determinant and making the change Zk = 
(zfe + t)/t, we hnally obtain: 

r,-x 



Z*°P.. ;ra = <fa'( N - 1) Yl trj 



3 = 1 



o 



() 



n 



1 



1) 



[(toj - Wk)(t 2 WjWk - 2AtWj + 1)] 



d s iu 
(27ri) s 



(3.6) 

Here Ci denotes a small, simple, closed, positively-oriented contour enclosing point 
z = l. 

Formulae (|3.1[) and (|3.6j) can also been derived by other methods (i.e., without 
using the equivalence of the algebraic and coordinate Bethe Ansatz), e.g., starting 
with vertical monodromy matrix formulation of Z^° p Ts , analogous to (|2.22j) for 
Z*!°\ rs , and next using the technique of paper [l8j to evaluate the matrix elementQ. 
We also mention that formula 



in a different form and for special values of t 
and A, has been found in the context of enumerative combinatorics IT51. 



4. Calculation of ^° t .. >rs 

Taking into account commutativity of B-operators, see (|2.9[) , and using relation 
(12.101) , we can obtain, in the usual spirit of the algebraic Bethe Ansatz calculation 
(for details see, e.g., [1]), the relation: 

A(X r ) f[ B(X P ) = £ TTlHrT II H X °» A ^ II B (^)MK). (4.1) 

13=1 a =l J{Aa ' Ar) f3=l 0=1 

Using this commutation relation and taking into account that 

N 

^+x,...,iv(A) \ltY +1 ,..., N ) = II «(A,^ fe ) |ty+x,...,iv> 5 (4-2) 

fe=s+l 



^We are indebted to P. Zinn-Justin for explaining us this alternative derivation. 



and also using (I2.13p , we obtain 



7bot 



E E 

a i — 1 CK2 — 1 



e n n *(^^)U 



a s =l j=l k=s+l 



n n /( a « 2 ,a^ 



= 1 fi^a^Kj) 

n /( a «-v) 



/3i=l 



/3 2 = 1 



j3 8 ^OL\ ,...,Ct s 



; z/ s+ i, . . . , v N ), (4.3) 

u N -s} ■= {l,...,AT}\{a 1 ,...,a s }. 



where {ai, 

To proceed further, it is convenient to introduce function 



IV (A) 



where d(X, A') := sin(A — A') and e(A, A') := sin(A — A' + 2rj). Expressing functions 
/(A, A') and g(X, A') appearing in (|4. 3[) in terms of functions ci(A, A') and e(A, A') and 
substituting the Izergin-Korepin determinant formula, see (I2.16[) . for the partition 
function standing in (|4.3I) . we arrive at the expression: 



(4.4) 



zbot 



lla=l Uk=a+1 a(\ a ,v k )b{\ ai v k ) 



Y\.i^<x<p^N d (A/3, A Q ) n s +i<j<fe<jv d{v 3 ,Vk) 



ai Ck.2 — 1 



<X s = l 
ct s ^ct\ . ..a s _ 1 



i 



det {^(Aa^i/g+fc)}. (4.5) 



, , e(A Q) , A Q J i<j,fc<Ar-s 

Here x(P> a ) = 1 if ;8 > a, and x(/?,«) = otherwise. 

Clearly, the multiple sum in (|4.5p reminds the Laplace expansion of some NxN 
determinant. This is also in agreement with the fact that since v r (X a ) = (a = 
r + 1, . . . ,N) all summations in (|4.5|) can be extended till the value N. To write 
down such a determinant formula, let us set A Q = A + £ a (a = 1, . . . , iV), where 
A is some arbitrary parameter, and parameters ,£/v are all different. Using 

again the fact that for a function /(£), regular near point £ = A, we can write 
/(A + £) = exp(£<9 e )/(A + e)| e= 0j we can bring (14. 5p to the form 



?bot 



ril^a<^^JV d ( X P> A «) IL+l<j<k«iV d i v h v k) 

exp(£i<9 £l ) ... exp(£i<9 e J </j(Ai,ivk) 
exp(£ 2 <9 £l ) ... exp(£ 2 <9 e J (/j(A 2 , f s +i) 



^(Ai,i/jv) 
<£>(A 2 , v N ) 



exp(^jv9 El ) ... exp(£;yd Ss ) ¥)(Ajv,i/ a+ i) 



x n ^ ( a ■ 



n 



i 



(4.6) 



ei = ...=e s =0 

We stress that this expression is valid for the inhomogeneous model. 
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Let us now perform the homogeneous limit. We regard A as the parameter of the 
weights of the homogeneous model, so that parameters £1, . . . , £jv and v s +i, ■ ■ ■ j V N 
are sent to zero in the limit. The procedure can be done along the lines of [3] and 
it is explained in full detail in As a result, we obtain the expression 

tp{\) i. 



rhat 



N{N-s) 



Sx— V(A) 1 
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2N-S-2 



(sinej)^ rj ' [sin(ej — 2r?)] r ^ 1 



[sin(ej + A - 7])] 



JV-s 



n 



V (A) C" 1 
1 



1 



sin(e,- ~e k + 2rf) 



...=E S =0 

(4.7) 

where in writing the determinant we have changed the order of columns, in com- 
parison with (|4.6[) . 

In order to represent (|4.7p in terms of a multiple integral, we first transform the 
N x N determinant representation (|4.7|) to an s x s one, given in terms of certain 
set of orthogonal polynomials. The construction is based on the following general 
facts from the theory of orthogonal polynomials (see, e.g., [3U] ). Let {P n ( x )}^=o 
be a set of orthogonal polynomials, 



P n (x)P m (x)fj,(x) dx = h n 8 n 



(4.8) 



where the integration domain is assumed over the real axis and we choose normali- 
sation such that P n (x) = x n + . . . , and let c n denote the nth moment of the weight 



x n /j,(x)dx (n = 0, 1, ...). 



(4.9) 



Then deti^- i fe^iv{cj+fe_2} = hohi ■ ■ ■ /i/v-i and, more generally, for s — 1, ...,N, 
the following formula is valid: 



CO 


Cl . 


■ CN-s-1 


1 


1 


1 


Cl 


Cl ■ 


CN-s 


X 1 


x 2 ■ 




Cjv-l 


C N . 


■ C2N-S-2 


x?- 1 


= h hi ■ 


■ • x s 

• hN-s-1 



det {Pjv- s +i-i(z fc )}. (4.10) 

In our case c n := d^ip(A), and the integration measure n(x)dx can be found through 
the Laplace transform for function ip(X); for explicit expressions, see |21) . 
As in [TT] . we denote 



where ip := </?(A), and ft.„ is as in 

a sine 

w(e) 



We also introduce functions 
b sine 



b sin(e — 2??) 
which, in particular, satisfy the relation 

b sin(e — 2?/) 



6(6) := 



a sin(e + 2rj) 



c sin(e + A — 77) w(e) 
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(4.11) 
(4.12) 
(4.13) 



where the short notations for the weights a :— sin(A + 77), b := sin(A — 77), and 
c := sin 27/ are used. Taking into account that ip — c/ab, and using the relation 

sin(ei + A + 77) sin(e 2 + A - 77) 1 (1 - u}(e 1 ))(u(e 2 ) - 1) 



sin(ei - e 2 + 277) (p J)(ei)w(e 2 ) - 1 

after applying (14. 10)) to (|4.7[) , we obtain: 

Z N /a\''i+---+ r 



(4.14) 



7bot 



s(2N-s + l) 3( 3 -3) 

a 2 2 c' 



(?) 



det {K N -,+j-i(d Ek )} 



><n 



n 



1 



a>(ej)w(£jfe) - 1 



=£3=0 
(4.15) 

In deriving of this formula, we have also used (I2.18[) to express a proper factor as 
the partition function Zn- 

Now we are ready to write representation (|4. 1 5[) as a multiple integral. We 
follow the procedure developed in |11) . The key relation here, valid for an arbitrary 
function f(z) regular at the origin, is 

1 r (z - I) 1 *' 1 
= 77- O ^ i h N {z)f(z) dz. (4.16) 

=0 27T1 J z" 

Co 

Here Co, as above, is a small, simple, closed, positively-oriented contour enclosing 
point z = 0, and h^iz) (not to be confused with h n in (|4.8p ) is the generating 

r( r ) „r-l 



K N ^(d e ) f(u(e)) 



function for the one-point boundary correlation function, hpf(z) = X^rLi H 



N 



where 



= K N ^{d e ) 



He) 



N- 



Me) - IP-* 



(4.17) 



e=0 



This function can be viewed as the s = 1 case of the row configuration probability, 

(r) (r) 

■ffjv := Hpj J v Indeed, in this case the partition function of the upper sublattice is 
simply Z' op = a N ~ r b r ~ 1 c, while Z^ ot can be found from (|4.15p . thus reproducing 

613). 

To write down the resulting multiple integral representation for Z^ ot Ts , we 
introduce functions 

deti^fc^gjzjT^fc ~ iy~ 1 h N -s+ 3 (z k )} 



hff, s (zi, ...,Z S ) 



(4.18) 



rii^j<fc^s( 2: fe z j) 

which can be viewed as multi- variable generalisations of hpf(z) (for a detailed dis- 
cussion of its properties, see [H]). Noticing that 

t 2 u{e) 



(4.19) 



2Atuj(s) - 1 ' 

where t = b/a and A = (a 2 + b 2 — c 2 )/2ab, we can readily rewrite the orthogonal 
polynomial representation (|4. 1 5[) in virtue of (|4.16p as follows: 



Z 



hot 



lS (N-l) c s 



II II 



Co Co 



■j = l J l^j<fe<s 



ZjZk — 2AtZj + 1 



x h NtS (zi,. . . ,z„) 



d"z 
{2iri) s 



(4.20) 
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This formula is the desired representation for 2y ot Ts , valid for the homogeneous 
model. 



5. Emptiness formation probability 



An important example of correlation function which can be built from the row 
configuration probability is the emptiness formation probability. As in we 
denote by Fj^' the probability of observing all arrows on the first s horizontal 
edges (counted, as usual, from the top) located between r-th and (r + l)-th vertical 
lines (counted, as usual, from the right) to be all pointing left. Equivalently, due 
to both the domain wall boundary conditions and the ice-rule, we can define it as 
the probability of observing the last N — r arrows between the sth and (s + l)th 
horizontal lines to be all pointing down, and hence (see also figure [5]) we have the 
relation: 

(5.1) 



7( r >«) 
JV 



E 



l^ri<r2<-"<r s ^r 

Our aim here is to address how this summation can be done for the row configuration 
probability, given by (|2.19[) . (|3.6[) and (|4.20[) . to reproduce the multiple integral 
representations for the emptiness formation probability obtained in 

Let us first recall the results of paper [TTj . The following two multiple integral 
representations have been obtained: 



7<>. s ) 
JV 



( ^ | •• J II 

C Co 1 = 1 
X flN,s(zi,- ■■,Z," 



[{t 2 - 2At)zj + 1] 



8 ~3 



d s z 



1>-J'+1 



n 



Zk 



t 2 ZjZk — 2 Ate,- + 1 



(27ri) s 



l a s(s-l) ( 



s'.a 



Co Co 



n 

3=1 



[{t 2 - 2At) Zj + l] s 



n 

j,k=l 



Zk 



t 2 ZjZk — 2Atzj + 1 



x h,N,s(zi,- ■ ■ , z s )h St s(u(zi), ■ • ■ , u(z s )) 



d s z 



where 



u(z) := — - 



1 



(5.2) 



(5.3) 



(t 2 -2At)z + l' 

The two representations in (|5.2[) are related by a symmetrization of the integrand, 
which is due to the following relation 



Asym 

2l,...,2 s 



-q [(t 2 - 2At)zj + l](t 2 ZjZ k - 2Aiz fe + 1) 



l^j<k^s 

z s 

,! a s(s-l) c s 



n 

3=1 



{Zj - 1) 
[(t 2 - 2At) Zj + l] 8 " 1 



Yl ( z k ~ Zj)h SiS (u(zi),...,u(z s )). 

l^.j<k^s 

(5.4) 



Ep(- 1 ) [PI /(^Pi J ---: z pJ, and the sum is 



Here Asym zu _ >Zt f(zi, . . . ,z s ) 
taken over permutations P : 1 , . . . , s i-> P\ , . . . , P s , with [P] denoting the parity 
of P. For details on the proof of relation (|5.4[) , see [TT] . 
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Let us discuss the representation for emptiness formation probability obtained 
from (|2.19[) . (|3.6I) and (|4.20[) . according to relation (|5.1[) . Direct substitution gives 

s r,— 1 



N 



O 



Ci 



( 



(27ri) ; 



o 



o 



n 



i 



Co Co ■ 



L («;,• - l) s 
1 v J 1 i<n 



e n 

<r2<--<r s ^r j=l 



>< n 

i<j<fc<« 



(wj — Wk){t 2 WjWk — 2Atwj + l)(zfe — zj) 



t 2 ZjZk — 2Atzj + 1 



J fiN, s (zi, 



,Z S ) 



d s z 



(2m) s 
(5.5) 

so performing here the multiple sum and integrating over a set of variables (e.g., 
over w%, . . . , w s ) we should reproduce, modulo symmetrization of the integrand, the 
s-fold integral representations (|5.2[) . 

To address this problem, let us first consider the evaluation of the multiple 
sum in (|5.5I) . Observing that the integral over Zj vanish for Tj 0, because in this 
case the integrand is regular at Zj = 0, we can replace the sum in (|5.5|) over values 
1 ^ n < r2 ■ ■ ■ < r s r with a sum over values — oo < r\ < T2 ■ ■ ■ < r s r. Then, 
denoting Xj — Zj/wj, the summation can be done using the identity 



s s 

E 11^ II 



1 



rr-s+j 



(5.6) 

\ - — \ ' I I I I ' \ , I 

-oo<ri<r2<---<r s sCrj — 1 

which can be easily verified by expanding the denominators in the right hand side 
in Taylor series. As a result, we find that (|5.5I) simplifies to expression: 



x — 



o 



Ci 



n . 



>< n 



(27Ti) s . 

Co Co 

(-Wj - w k )(t 2 WjW k - 2Atwj + l)(z fc - zj) 
t 2 ZjZ k — 2Atzj + 1 



h N ^(zi,...,z s ) 



(2?ri) s 
(5.7) 



and we are left with performing an s-fold integration. 

We shall integrate over variables W\ , . . . , w s in (|5.7[) . Let us consider the equiv- 
alent integral where the integrand is symmetrized with respect to permutations of 
these variables. Define function 



,w s ;zi, 



■ z s ) 



n 



i 



Wk - Wj 



x Asym 

Wl,...,1D s 



(5. 



Integration over the Wj 's is done with the result 



n 



7 T~T 9 

j=i v J ' isjj<fc^s 



d s w 
(2ni) s 



= (-l)"(- 1 )/ 2 «! #.(!,..., 1;^,...,%), (5-9) 



which can be easily found by noticing that in evaluating the residues one has to 
differentiate only the factor Ylj <k ( w j — w k) 2 ■ 
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Finally, the desired result for the emptiness formation probability amounts in 
proving the identity: 



s! Asym 

Zl,...,2 s 



2Mz k + 1)] 

r±y(s+i)/2 



Asym 



n 



[(t 2 - 2M)zj + l}(t 2 ZjZ k - 2Mz k + 1) 



(Zj - 1) 



(5.10) 



This identity has to be used together with identity (|5.4[) to reproduce (j5.2|) . We 
find identity (|5.10p rather difficult to prove directly, and presently we have only 
been able to verify it through computer-aided calculations for small values of s. We 
note that rather similar identities have been discussed in 

In conclusion, in this paper we have introduced and calculated a nonlocal cor- 
relation function of the six-vertex model with domain wall boundary conditions, 
the row configuration probability. It is given as a product of two factors which can 
be treated as the partition functions on upper and lower sublattices of the original 
lattice (see figure [5]). We have represented these partition functions in terms of mul- 
tiple integrals, see (|3.6|) and (|4.20l) . The row configuration probability can be used 
for computing other correlation functions, provided that sums like those appearing 
in (|5.1j) can be evaluated. To illustrate this, we have considered the problem of 
reproducing the known result for the emptiness formation probability. We have 
shown that in this case the problem boils down to identity (15.101) . A direct proof 
of this identity, in addition to the indirect one following from the known equality 
of (|5.2p and (|5.5p . could be useful for the evaluation of other correlation functions. 
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